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1 INTRODUCTION 



ABSTRACT 

We investigate the potential of weak gravitational lensing maps to differentiate be- 
tween distinct cosmological models, considering cosmic variance due to a limited map 
■ extension and the presence of noise. We introduce a measure of the differentiation 

between two models under a chosen lensing statistics. That enables one to determine 
in which circumstances (map size and noise level) , and for which lensing measures two 
models can be differentiated at a desired confidence level. As an application, we sim- 
Cn| \ ulate convergence maps for three cosmological models (SCDM, OCDM, and ACDM), 

calculate several lensing analyses for them, and compute the differentiation between 
the models under these analyses. We use first, second, and higher order statistics, 
including Minkowski functionals, which provide a complete morphological character- 
ization of the lensing maps. We examine for each lensing measure used how noise 
affects its description of the convergence, and how this affects its ability to differenti- 
ate between cosmological models. Our results corroborate to the valuable use of weak 
^ ■ gravitational lensing as a cosmological tool. 

\ Key words: gravitational lensing — large-scale structure of universe 

6 

, The mass inhomogeneities in the Universe leave an imprint in the light traveling through it: that is gravitational lensing. The 

retrieval of this information can be of singular value to narrow down cosmological models. 

The statistical shape distortion of distant galaxies, or cosmic shear, is is one of the aspects of gravitational lensing in the 
weak regime (small deflection angles), and has been observed by several groups (see Mellier 1999 for a review up to this year, 
and for more recent measurements see Van Waerbeke et al. 2000; Wittman et al. 2000; Fischer et al. 2000; Bacon, Refregier 
& Ellis 2000; Kaiser, Wilson & Luppino 2000; Maoli et al. 2001; Rhodes, Refregier & Groth 2001; Van Waerbeke et al. 2001; 
Hoekstra, Yee & Gladders 2001). Concomitantly, theoretical works in the area indicated that such measurements can be very 
revealing about the large-scale structure of the universe (Kaiser & Squires 1993; Bernardeau, Van Waerbeke & Mellier 1997; 
Jain, Seljak & White 2000; Bartelmann & Schneider 1999, 2001), in addition to the study of individual galaxy clusters. The 
main attractive of gravitational lensing being that it probes mass directly, avoiding issues such as mass-light bias. 

The perspective of a growing number and quality of weak lensing measurements stimulates the question of how far the 
use of weak gravitational lensing as an astrophysical and cosmological tool can go, and how much can be achieved. That 
ultimately depends on the extension and quality of the lensing maps, the amount of information that can be extracted from 
them, and the ability of this information to differentiate between theoretical models. 

The aim of this paper is to investigate the potential of weak gravitational lensing maps to differentiate between different 
cosmological models. We incorporate limitations resulting from the presence of noise and cosmic variance for several lensing 
measures, including a morphological analysis, the Minkowski functionals (Matsubara & Jain 2001; Sato et al. 2001). 

Theoretical predictions for lensing measures can be obtained by two approaches. The first is the analytical, which is based 
on gravitational lens theory, and resorts frequently on perturbation theory to calculate expressions for various statistics in an 
assumed model. The second approach is to simulate realizations of lensing maps in a chosen model, and directly measure the 
quantities of choice. Here we review the first, and use the second for an application example. 

We obtain that even for modest field sizes (9 degrees 2 ), and noise at the level of current surveys, the cluster normalized 
cold dark matter models (SCDM, OCDM, and ACDM) can be differentiated with significant confidence. Second order statistics 
(convergence variance and lensing measures dependent on it) proved to be the best discriminatory analyses for these models. 
However, we also obtain that higher order statistics (measures that are sensitive to non-Gaussian features of lensing maps) can 
differentiate between some models in certain circumstances. This result is particularly important when aiming to differentiate 
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cosmological models that have the same power spectrum, but distinct non-Gaussian properties. We suggest that Minkowski 
functionals should be included in the row of statistics normally used to extract information from lensing maps, because 
their morphological characterization of maps contains information that is independent of measures such as the probability 
distribution function of the convergence. 

The paper is organized as the following. In Section £j we briefly review how lensing maps can be generated and analyzed, 
and we introduce a measure of the differentiation between two models through a chosen lensing statistic. In Section w we 
compare three cosmological models (SCDM, OCDM, and ACDM) as saw through weak gravitational lensing. Our conclusions 
are in Section m. 



2 LENSING MAP 
2.1 Generation 

The image distortion of background galaxies by weak gravitational lensing can be expressed by the distortion matrix (Bartel- 
mann & Schneider 2001) 

A(9)=( 1 -*- 71 -\ ), (1) 

V -72 1 - K + 71 J 

where k is called the convergence, and 7 = 71 + 172 = \^f\e 2l,p is the shear. 

Matrix A transforms a circular source image in an ellipse with semi-axis stretched by a factor (1 — n ± |7|)~ from the 
original radius, and magnification /j, — (det.4) -1 = [(1 — k) 2 — 7 2 )] -1 . 

The shear can be directly estimated from ellipticity measurements of background galaxies, and shear and convergence 
can be mapped in each other through 



7(0) = ~ j w(# - e') K (e')d 2 e' , 



(2) 



with 



or through the relation between their Fourier transforms 

ii — 12 ~ m . 2Z1Z2 . 

p-7i(0 + -p ' 

We can compute the components Aij of the distortion matrix by calculating the angular deflection of a photon traveling 
from a comoving distance w through a gravitational potential <E> 



*(9 = ^V^W + n^T*(0 • (4) 



Aij(0) = -2 g(w, w)d l d J ;$(0, w')dw , (5) 



where 

/ / \ fK(w')f K {w - w') 

9{w > w) = Mw) ' (6) 

and /k(w) is the curvature-dependent comoving angular diameter distance. 
From equation (||) and using Poisson equation, 

V^=l H lnJ-, (7) 
where 5 is the density contrast and a is the scale factor, we can write the convergence as 

K{0) = f g( w \ w ) 5 J^O-dw' . (8) 

2 J Q a(w>) 

The convergence is a weighted projection of the mass inhomogeneities between source and observer. 

In this work we concentrate on the convergence field k(0) as our lensing map. It can be derived from a real or simulated 
shear map through equation (^), or computed directly from a N-body simulation of 8(0, w) and the use of equation (^). We 
choose the second approach for our application example at Section ^. 



2.2 Statistics 

Here we review some statistics used to characterize the convergence field, the angular power spectrum, second and higher 
order statistics, the probability distribution function (PDF), and Minkowski functionals. 

The convergence field can be expanded in Fourier modes, and their amplitude averaged for each wavelength. That is the 
direct calculation of the angular power spectrum of the convergence from a real or simulated lensing map, 
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p«{1) = (l«(OI 2 ) 



(9) 



An analytical prediction for the convergence angular power spectrum can be obtained from equation (g) , and be expressed as 
a weighted integral of the time-evolving density power spectrum Ps (k, w) 



2nl 2 P K (l) = 36TT 2 n 2 m l 2 



g 2 (w,w ) 



\w)P s 



k = 



I 



dw 



(10) 



However, the angular power spectrum does not contain all the statistical information about the Ac-map. All the information in 
the phases of the complex Fourier modes is lost. These phases carry the non-Gaussian features of the map, which are important 
signatures of non-linear evolution, and also of models of structure formation containing non-Gaussian initial conditions (as 
opposed to most of the inflationary scenarios). 

The Fourier transform of the angular power spectrum of the convergence is the angular two-point correlation function 
Cre(r) = (k(6)k(0 + r)). The convergence field variance a 2 = {k 2 ) can be seen as the value of the angular two-point correlation 
function at the origin, a 2 — C K (0), and expressed as (Jain & Seljak 1997; Bernardeau, Van Waerbeke & Mellier 1997) 



(k 2 ) = 36tt 2 Q. 2 



kdk 



g 2 M 

a 2 (w) 



■Ps{k,w)Wi[kf K (w)6 s ]dw , 



(11) 



where W2 = 2Ji(x)/x is the Fourier transform of the top-hat window function (Ji is the Bessel function of first order). 
Expressions for higher order powers, (k 3 ) and (k 4 ), can be similarly obtained as integrals of the bispectrum and trispectrum, 
respectively. See Hui (1999), and Munshi & Jain (2000), or Cooray & Hu (2001) for more elegant expressions. Related statistics 
are the skewness S = (k 3 )/o" 3 , the kurtosis K = ((k 4 )/(j 4 ) — 3, or the more general moments Sn = (k N )/<t« . 

The probability distribution function contains in principle more information than the variance and higher order moments 
of the convergence field taken individually, which can be calculated from the PDF itself, 



(k n ) = I F(k)k ' dn . 



The PDF can be expanded around a Gaussian through the Edgeworth approximation (Jusziewicz el al. 1995) 



F(k) = 



2 /2a 



l + a K ^R 3 ^/a K ) + a 2 K 
2tv(t k I ol 



|iH 4 ( K /a K ) + | 3 -H 6 ( K /a K 



+ 



(12) 



(13) 



where H,, is the Hermite polynomial of order r\. 

Minkowski functionals contain all the morphological information about a convex body, and for a Gaussian field they can 
be calculated exactly (Winitzki & Kosowsky 1998). That makes the Minkowski functionals of the convergence maps a very 
interesting statistic to use (Matsubara & Jain 2001; Sato et al. 2001). A threshold value v defines excursion sets in which the 
value of v(&) = n(G)/a K is larger than the threshold: vo(v) is the fractional area of the map above the threshold, v\(v) the 
boundary length (per area), and V2(y) the Euler characteristic. This last functional is equivalent to the topological genus of 
the map, roughly speaking, the number of disconnected regions minus the number of holes. 

Approximate expressions for the Minkowski functionals can be obtained through perturbation theory about the exact 
expressions for a Gaussian field (Matsubara 2000) 



V2{v) 



2 \V2~J 6 v 7 ^ 
-e ' i 1 + a K 



8^^ 



,«)) 



-,(D 



-|~H 3 (^) + -|-Hi(i/) 



1 <f\ -u 2 /2 I,,, 

— e ' { Hi(i/) + a K 



2(2tt)t ol 



,(()) 



Oct 1 ) Q 1 



3 



(14) 



(15) 



(16) 



where a\ = ((Vk) 2 ), and skewness parameters are defined by S^ u> = (k 3 )/ct 4 , 
-3((Vk- \7n)(\7 2 K,)}/at. 



-(3/4:)(K 2 (V 2 K))/(ala 2 ) and Sf > 



2.3 Differentiation 

Here we introduce a measure aimed to quantify the difference between two models, or a model and an observational result, 
under a chosen lensing map analysis (Guimaraes 2001). This quantity is also aimed to allow the comparison of different lensing 
statistics, for instance, being able to point which analysis is the most appropriate to use when trying to differentiate between 
models through weak lensing. We construct one such quantity in Appendix [AL 

The differentiation T> between models A and B under the lensing analysis Y is defined as 

V[Y] = 1 - e-* 2/2 , (17) 
where 
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Figure 1. Convergence maps for SCDM (left), ACDM (center), and OCDM (right). The top maps are pure (no noise added), and the 
bottom maps have noise added. All maps are 3x3 degrees 2 , and smoothed at 1 arcmin. 



Table 1. Cosmological models considered. Other parameters: h = 0.7, Y = Q m h, Hf, = 0, L^ ox = 128h 1 Mpc, l ao f t = 30h 1 kpc, 

Npart = 128 3 . 





" 7 ? J 


n A 




SCDM 


1.0 





0.56 


OCDM 


0.3 





0.84 


ACDM 


0.3 


0.7 


0.99 



v 2_ 1 ^ [Ya(p>)-Y b ( Pi )] 2 

x - n2^ G 2 A(Pi) + al{pi) ■ w 

Ya(Pi) is the mean value of the lensing measure Y for model A at N values of the analysis parameter p, an d its variance is 



a A (pi). Equivalently for model B. Equation ( |18[ ) is a particular discrete form of the more general expression (All), when the 
parameter interval is divided in N equal-size segments. 

According to this definition, the differentiation is quantified on a scale from to 1. For models that are similar under a 
given analysis T>[Y] w 0, and for very distinct models T>[Y] ~ 1. If one wants to study a set of models, it is desirable to find 
which lensing measures give the largest differentiations, so the discrimination between models is robust. On the other hand, 
if one has an observational result, the interest is to find which model gives the lowest differentiation for the chosen analysis. 
That is the best fit model. 

If we assume that the differentiation has a Gaussian distribution, then it follows that its variance is 

A 2 [V] = X V* 2 . (19) 



3 STUDY OF WEAK LENSING MAPS 
3.1 Map construction 

We compare three important cosmological models, SCDM, OCDM, and ACDM (see table]]]). We used the Hydra N-body code 
(Couchman, Thomas & Pearce, 1995) to simulate the large-scale mass distribution in these models by evolving 128 3 particles 
from a redshift z = 50 to z — 0. The particle coordinates at juxtaposed boxes from a redshift z — 1 (our assumed source 
redshift) to z — were saved, and projected into a middle plane in each box (the orientation and origin of the projections were 
randomized). The resulting set of 2D density contrast fields (in 1024 2 grids) at those positions in redshift space were used as 
a discrete approximation for the computation of the convergence field according to equation ([s]), that is usually refered to as 
multiple-plane lens method. 

The main source of noise in real lensing maps is the intrinsic ellipticities of the background galaxies. We simulated this 
noise by adding a Gaussian random field to the pure convergence following Van Waerbeke (2000), which showed that this is 
a good approximation. We then used a top-hat window of radius 8 S to smooth the map, which can be described by 

k(0) =K o (0) + n{0) , (20) 

where Ko(0) is the smoothed pure convergence, and n(9) is the smoothed noise field. The noise part is a Gaussian correlated 
field of variance cr 2 , 

(21) 



2n g n9'i 



0.16 was the intrinsic ellipticity variance adopted, and n a — 60 arcmin - was the mean source galaxy density 
assumed. The top-hat smoothing of the original Gaussian random field introduces correlations in the noise at scales bellow 
twice the smoothing radius. 

Figure |l| shows some realizations of the K>maps. We generated 25 realizations for each model, however, because we only 
used one N-body simulation for each model, the K-maps generated cannot be considered to be totally independent. That is 
not a source of major concern for our objectives in this work. 
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Figure 2. Top panel: angular power spectrum of the convergence field smoothed at 0.25 arcmin scale. Bottom panel: angular power 
differentiation between models ACDM-SCDM (dashed lines), and ACDM-OCDM (dotted lines). Thick lines are for pure maps, and thin 
lines for noisy maps. Only the error-bar tips are show for clarity. 



3.2 Statistics and differentiation results 

We applied a number of statistical measures to the convergence maps of the three cosmological models. We considered pure 
(no noise) and noisy maps of field sizes 3x3 degrees 2 , and lxl degree 2 . For each analysis we calculated the differentiation T> 
according to equation @ between ACDM and SCDM, and between ACDM and OCDM. 

Figure ^ shows the results for the angular power spectrum of the convergence for a field smoothed with 6 S = 0.25 arcmin, 
and its differentiation between models. The total power of the convergence field, P K , is just given by the sum of the power of 
the pure field P Ko and the power of the noise field P n (a power law), because the two fields are uncorrelated (by construction), 

P K {l)=P Ko (l)+P n (l), (22) 

and its variance can also be fragmented as 

A 2 [P K ] = A 2 [P Ko ] + A 2 [P n ] + A 2 [<K n*)] + A 2 [<^n)] . (23) 

For low values of the wavenumber I (large scales) noise gives a small contribution to the power of the convergence, but the 
variance of the measured power is large due to the restricted sampling. This large variance implies a lower differentiation 
between models at low I. For high I (small scales) the measured power is suppressed by the field smoothing at wavenumber 
values above l s ~ d^ 1 . The noise field gives a major contribution to the total power and power variance, such that the 
examined models become indistinguishable as quantified by the low values of T> [2nl 2 P(l)] in this case. Note that the analyzed 
maps have different areas (ACDM 9.6 degrees , SCDM 15.6 degrees 2 , and OCDM 12.3 degrees 2 ), therefore the power variances 
cannot be directly compared. 

The addition of noise transforms the probability distribution function for the pure convergence field F Ko according to the 
convolution 

F K {x)= / F Ko (y)^= dy. (24) 



For a small smoothing angle 9 S the noise variance cr 2 is large, so the convolution has a large effect on the pure convergence 
PDF. fn contrast, for a large 8 S , cr 2 is small, and F K does not differ substantially from F Ko . This effect is illustrated on Figure 
^ which shows the results for the probability distribution function of the k field, and the integrated PDF differentiation 
in the interval —0.04 < k < 0.04. The differentiation curves indicate that noise reduces the differentiation between models, 
but moderately. Smoothing has a much more visible role in reducing this differentiation - the convergence maps for different 
models are made more homogeneous by smoothing, and therefore made more alike. The map size also has a major effect on 
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Figure 3. Top panels: convergence probability distribution function for pure (left) and noisy (right) maps at 1 and 10 arcmin smoothing 
scales of 3x3 degrees 2 fields. Bottom panels: PDF differentiation between ACDM-SCDM (bottom left) and ACDM-OCDM (bottom 
right), as function of the smoothing angle 8 S . Thick lines are for pure maps, and thin lines for noisy maps; solid lines are for 3x3 degrees 2 
fields, and dotted lines arc for lxl degree 2 fields (the superior error-bars for these are not show for clarity). 



the ability of the PDF analysis to differentiate between models, because the variance of the PDF measurement increases with 
a reduced field size. 

Figure ^ shows the convergence variance (k 2 ) as a function of the smoothing angle, and its model differentiation. Because 
the pure convergence and noise fields are uncorrelated, the convergence variance can be written as 

< K 2 ) = ( K 2 ) + al , (25) 

A> 2 ] = A 2 [a 2 Ko ] + A 2 [al] + 2A 2 [{k u}] . (26) 

Noise dominates for small field smoothing, because a 2 , oc 0~ 2 , while the variance of the pure field decreases with a smaller 
power (cr 2 ~ 97°' S ), and is negligible for large 6 S (error-bars are increased though). The curves for T> [(k 2 )] demonstrates 
that the field variance is a powerful analysis to discriminate between the models considered, even for small maps. 

The mean third power of the field is a direct measure of its deviation from Gaussianity, and its value is not affected by 
the noise field (see Figure ^), because (n) — (n 3 ) = 0, 

<* 3 > = <*o> . (27) 
Unfortunately the noise field increases (k 3 ) error-bars, reducing the differentiation between models 

A 2 [( K 3 )] = A 2 [( K 3 )] + A 2 [(n 3 )] + 2A 2 [( K 2 n)] + 2A 2 [< Ko n 2 )] . (28) 

Figure ^ also shows the results for the related quantity skewness, S = («; 3 )/(T 3 . Its null differentiation between ACDM and 
OCDM, and not null between ACDM and SCDM, suggests that the convergence skewness is sensitive to the mass density, 
and insensitive to the cosmological constant. 

Figure shows the results for the mean fourth power of the convergence field (k 4 ), and kurtosis K = ((n 4 }/at) - 3. The 
contribution of the noise field can be visualized in the expression 

(K 4 ) = (4} + (n}+6(K 2 n 2 ) , (29) 

where (n 4 ) = 3<r 4 . This coupling of the pure field with the noise field, {K 2 n 2 } term, allows the model differentiation for (k 4 ) 
to be larger for the noisy maps than for the pure ones. The kurtosis of Gaussian fields is null, and that is what is observed for 
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0.1 1.0 10.0 0.1 1.0 10.0 

8 S (arcmin) 8 S (arcmin) 

Figure 4. Convergence variance of 3x3 degrees 2 fields (top panel), and convergence variance differentiation between ACDM-SCDM 
(bottom left) and ACDM-OCDM (bottom right), as function of the smoothing angle 8 S . Thick lines are for pure maps, and thin lines 
for noisy maps. In the bottom panels solid lines are for 3x3 degrees 2 fields, and dotted lines are for lxl degree 2 fields (the superior 
error-bars for these are not show for clarity). 

noisy convergence maps at small and large smoothing angles. At small 8 S because noise dominates, and at large 6 S because of 
the extreme field smoothing. 

Some Minkowski functionals curves for the convergence field are shown at Figure [?], and the differentiation results for 
this analysis of the models are shown at Figure || Minkowski functionals are very sensitive to noise and smoothing. The 
morphology of the convergence maps for the examined models seem to not differ greatly as described by the functionals. 
However for noisy fields the second and third functionals have considerable discriminatory power at small smoothing angles. 
It is in principle counter-intuitive that noisy maps would have a larger differentiation than pure convergence maps, but as 
seen in the case of the (k 4 ) analysis, noise can introduce a "statistical contaminant". The discussion on Minkowski functionals 
of noisy maps presented in Appendix ^ may be useful in understanding this effect. 

Another interesting result to note is the low differentiation between models under the first Minkowski functional analysis 
(fractional area), which can be saw as a cumulative PDF. Figure ^| shows a strong differentiation between models under 
the PDF analysis, so the results for T>[vq] could seem contradictory with those for T>[PDF], however when we remember 
that Minkowski functionals are displayed as functions of the threshold v — k/ct k this apparent discrepancy disappears. This 
observation suggests that the large PDF differentiation between the models is due primarily to their very distinct convergence 
variances. 



4 CONCLUSIONS 

We study the potential of weak gravitational lensing maps to differentiate between different cosmological models, using a 
variety of statistical measures, including morphological analysis, and taking in account the effects of noise and the uncertainties 
resulting from a limited field size. 

We reviewed how gravitational lensing maps can be generated, and also revisited some statistics that can be used 
to characterize the lensing maps, or extract information from them. We introduced the quantity T>, which quantifies the 
differentiation between two sets of lensing maps under a given statistic. We used simulations of convergence maps in three 
cosmological models (ACDM, SCDM, and OCDM) to study and compare these models. Several of the analyses considered 
here were also investigated in previous works, with results consistent with ours. The novelties presented in this paper are the 
quantitative description of the differentiation between models (which proved to be very useful), the possibility to compare 
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Figure 7. Minkowski functionals for the convergence field (pure and noisy), for smoothing angles 0.25 (left) and 4 (right) arcmin. The 
first functional vq is dimensionless, v\ units are arcmin" 1 , and vi arcmin -2 . 
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Figure 8. Minkowski functional differentiation between ACDM-SCDM (left) and ACDM-OCDM (right), as function of the smoothing 
angle 8 S . Thick lines are for pure maps, and thin lines for noisy maps; solid lines are for 3x3 degrees 2 fields, and dotted lines for lxl 
degree 2 fields (the superior error- bars for these are not show for clarity). 



side-by-side different statistics, and a systematic consideration of the role of noise and field size in all lensing measures. It was 
also the first time that the Minkowski functionals for noisy maps were calculated. 

Our results show, as expected, that the lensing measures of small fields have large variance, because of limited sampling. 
These large variances imply a small value for the differentiation T>: the observation of a too small sky patch does not allow 
the discrimination of cosmological models. Other general behavior, independent of the analysis or model considered, is the 
approximation of the differentiation values obtained for pure and noisy convergence maps when the smoothing angle is large. 
As smoothing increases the noise variance decreases (noise becomes irrelevant for a sufficiently high smoothing angle), but 
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the pure convergence is also homogenized, and the differentiation between models becomes small. Our results suggest that it 
is advantageous for model discrimination purposes to use a minimum smoothing, even considering that this implies a noisier 
map. 

For purposes of discriminating the three models examined in our simulations, the variance and the angular power spectrum 
of the convergence (two of the most simple and popular analysis) proved to be very good lensing measures. Even for small 
field sizes, and noisy maps, the convergence variance can differentiate between the models at a great confidence level (if the 
smoothing is not too extreme). 

Higher order statistics ((fi: 3 ), skewness, (ft; 4 ), kurtosis, and Minkowski functionals) have higher variances (error bars) in 
general, and lower discriminating power for the examined models. However, the information obtained from these statistics is 
precious for its independence from the information offered by statistics such as field variance, and angular power spectrum. It 
is possible to imagine two cosmological models that have the same mass power spectrum, but different non-Gaussian features. 
We found that the second Minkowski functional (vi) is an equal or better model discriminant than the more commonly used 
third functional (the topological genus), and is very competitive in relation to other measures. This suggests that Minkowski 
functionals should be included in the row of available statistics to maximally extract information from weak lensing maps. 
They should be particularly useful to differentiate models that have the same PDF, but distinct morphology. 

The presence of noise makes the extraction of information about non-Gaussian features of the pure convergence much 
more intricate. While for the analyses that are blind to non-Gaussianities (variance and angular power spectrum) the noise 
term is clearly isolated from the pure convergence term, for analyses sensitive to non-Gaussianities (wit h exception of (k 3 )) 
the noise terms are entangled with the pure convergence terms (see equations ^3, ^jj], [B2|, and [B5-B7]). So, even for a 
simple noise model as the one used in our simulations, the modification of the lensing measures due to noise is nontrivial. That 
implies that the retrieval of information about the pure convergence field (and ultimately, from the underlying large-scale 
structure and cosmological model) from noisy maps is also nontrivial - even more when the noise field is not well known. 
Minkowski functionals are a particularly severe example of this entanglement between the pure convergence with noise (as 
shown in Appendix [b]) . 

The complicated analytical description of lensing measures sensitive to non-Gaussian features of noisy maps, and the 
intrinsic difficulty of incorporating observational aspects trough the analytical approach, favors the use of simulations over 
an analytical approach for observational predictions, or parameter estimation from real lensing maps. These arguments add 
to the work of Jain, Seljak & White (2000), which points to the limitations of perturbation theory in providing accurate 
predictions for most weak lensing statistics at small scales (low smoothing). 

Different statistics of weak gravitational lensing maps reveal distinct features of its originating cosmic mass distribution 
and geometry. Therefore, a comprehensive and realistic study of the underlying cosmological model through weak lensing 
requires the use of a set of statistics, and the understanding of how these statistics and their variances are affected by 
observational constraints such as field size and the presence of noise. 
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APPENDIX A: CONSTRUCTION OF THE DIFFERENTIATION 
Al Simple case 

Let A and B be models or observations, Y an analysis of lensing maps, and a 2 its variance. Ya and Yb are the result of 
applying the analysis Y on maps of A and B. Ya is assumed to follow a normal distribution of mean value Ya and variance 

F[ Y A ] = —^ e ~ {YA ^ A)2/2 ^. (Al) 

Equivalently for B. 

We define the analysis difference between the two models as 

D = (A2) 
It is easy to see that D is a unit normal distribution, 

F[D] = -L e - {D ' D)2/2 , (A3) 



where the expectation value of D is 

5=4=% (A4) 

\J a + °i 

In fact, we are interested in the absolute value of D, so the distribution function for \D\, is 
F[\D\] = F[D]+F[-D] (A5) 

-(D-D) 2 /2 + e -(D+D) 2 /2 



If D — we say the two models are similar. So the probability of D — gives a good measure of how similar the two 
models are under the analysis Y: 

^l(0) = yfe- B2/2 , (A6) 

which is normalized to 1 over integration in D (from to oo). 

Two very different models (D 3> 0) have F\ D ^(0) w 0; in contrast, two similar models have a large probability that 
\D\ — 0. The maximum possible value for Fmi(0) is obtained when D — 0, which implies 

max[F p |(0)] = J| (A7) 

We can finally define the differentiation between the two models under the analysis Y as 

V[Y] = 1- F \ Dli °] =l~e- 52 / 2 . (AS) 
1 J max[F| D |(0)] V ; 

A2 Generalization 

We can generalize the differentiation for statistics that are functions of a parameter p (e.g. the angular power spectrum as 
a function of the wavenumber !, or a Minkowski functional as a function of the threshold v). Note that p, and the analysis 
itself, can be a vector. 

The same elements used before for the simple case can be generalized to Ya,b(p), ca,b(p), 
%) = ^ft^ (A9) 
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Using a reasoning similar to the one used in the simple case, we can obtain 

V[Y] = 1 - e~ x2/2 , (A10) 
where 



2 

X 



W(p) dp 



D 2 (p)W(p)dp, (All) 



and W(p) is a window function (which limits the parameter domain). 

Note that the discrete form of \ 2 assumes a very recognizable representation (equation j^]). It becomes then clear that 
the maximization of T>[Y] is equivalent to the minimization of x 2 , best know as the least square method. 

The differentiation can be further generalized to include the effects of covariance by 



v 2 

K N 



1\ 

V 2[$U(pi) - YBip^M-'MYAiPj) - Y B { Pj )] , (A12) 



where 



My = ([Y A (p t ) - Yb (pi)} [Ya( Pj ) - Y B (p 3 )}) (A13) 

is the covariance matrix, and the brackets denote ensemble average. Note however that we did not include the effects of 
covariance in our calculations in Section H. 



APPENDIX B: MINKOWSKI FUNCTION ALS OF NOISY MAPS 



The problem of expressing the Minkowski functionals for the noisy convergence map k in terms of the functionals for the pure 
convergence «o and the smoothed noise n seems to be an open issue. For the first functional, the fractional area, we found a 
solution, which can be obtained when one makes use of the fact that this functional is the cumulative probability function 



v (v) = 



F(x)dx , 



(Bl) 



and writes the convergence probability distribution function F K as the convolution of F K0 and F„ (equation pi| ), obtaining 



W) 



(n) 



J K o) 



dv . 



(B2) 



Another approach is to use the approximations ( |l^ - |l^ ) and calculate how their terms change for a noisy map. In fact, 
these approximations should be better representations of the noisy map than of the pure map, because noise addition makes 
the convergence closer to a Gaussian field. 

The term a{ = ((Vk) 2 } can be decomposed as 

<(V K ) 2 ) = <(V^ ) 2 ) + <(Vn) 2 ), 



and it can be useful to note that 

(Vk(0)v«(0')) = -c"(\e' - e\) , 



where C"(r) is the second derivative of the two-point correlation function. From (B4) follows that ((Vre) 2 } = — C"(0). 
The skewness parameters S^ ' is altered by the addition of noise by 

c(o) _ (^ 3 ) 
3 (a 2 + oiy ' 



2^V72 \\ / 2/V72 \\ i / 2/v72„\\ i / 2/T-72 . t-,2 



and to obtain S3 and S3 one needs to calculate 



(k 2 (V 2 k)) = {^(V 2 k )) + <^(V 2 n)) + {n 2 (V 2 K + V 2 n)) , 
((Vk • Vk)(V 2 k)) = ((Vko ■ V«: )(V 2 ko)) + ((2V«o ■ Vn + Vn ■ Vti)(V 2 k )) + ((V« + Vra) • (V« + Vn)(V 2 n)) 



(B3) 
(B4) 

(B5) 

(B6) 
(B7) 
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